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The rate of gross biological dimethylsulfide (DMS) production at two coastal sites west of the Antarctic 
Peninsula, off Anvers Island, near Palmer Station, was estimated using a diagnostic approach that 
combined field measurements from 1 January 2006 through 1 March 2006 and a one-dimensional 
physical model of ocean mixing. The average DMS production rate in the upper water column (0-60 m) 
was estimated to be 3.1 + 0.6 nM d 1 at station B (closer to shore) and 2.7 + 0.6 nM d" 1 at station E 
(further from shore). The estimated DMS replacement time was on the order of 1 d at both stations. 
DMS production was greater in the mixed layer than it was below the mixed layer. The average 
DMS production normalized to chlorophyll was 0.5 + 0.1 (nMd ^/(mgm -3 ) at station B and 
0.7 + 0.2 (nM d ')/(mg m ~ 3 ) at station E. When the diagnosed production rates were normalized to 
the observed concentrations of total dimethylsulfoniopropionate (DMSPt, the biogenic precursor of 
DMS), we found a remarkable similarity between our estimates at stations B and E (0.06 + 0.02 and 
0.04 + 0.01 (nM DMS d~')/(nM DMSP), respectively) and the results obtained in a previous study from 
a contrasting biogeochemical environment in the North Atlantic subtropical gyre (0.047 + 0.006 and 
0.087 + 0.014 (nM DMS d '^/(nM DMSP) in a cyclonic and anticyclonic eddy, respectively). We propose 
that gross biological DMS production normalized to DMSPt might be relatively independent of the 
biogeochemical environment, and place our average estimate at 0.06 + 0.01 (nM DMS d ’)/(nM DMSPt). 
The significance of this finding is that it can provide a means to use DMSPt measurements to extrapolate 
gross biological DMS production, which is extremely difficult to measure experimentally under realistic 
in situ conditions. 

© 201 1 Elsevier Ltd. All rights reserved. 


1. Introduction 

Dimethylsulfide (DMS) is a climate-relevant gas produced in 
the marine food web from the phytoplankton-derived precursor 
dimethylsulfoniopropionate (DMSP). Compelling observational 
and modeling evidence suggests that oceanic DMS emissions are 
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the largest biogenic source of sulfur to the atmosphere (Kettle and 
Andreae, 2000; Lana et al., 2011). DMS-derived aerosols influence 
the global radiation budget directly by scattering incoming 
photons and indirectly by acting as cloud condensation nuclei, 
thus increasing cloud albedo (Charlson et al., 1987; Ayers and 
Gillett, 2000; Leek and Bigg, 2010). DMS-derived aerosols are 
particularly important in the Antarctic region, due to high DMS 
concentrations in the surface waters during the austral spring and 
summer, coupled with low aerosol concentrations in the Antarctic 
troposphere, remote from terrestrial aerosol sources (Vallina 
et al., 2006). DMS is highly reactive in the atmosphere and has 
a short lifetime; consequently, atmospheric DMS concentra- 
tions are strongly influenced by local DMS emissions from the 
ocean. Ocean-atmosphere DMS fluxes, in turn, depend on the 
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sea-surface DMS concentrations, which are controlled by biological 
processes involving the entire marine food web, as well as physical 
and chemical processes in the water column. One of the central 
problems in quantifying and modeling the DMS-climate connection 
is the lack of understanding of the dynamics of the marine DMS 
cycle in the upper water column (Simo, 2001). In particular, the 
most fundamental quantity in the budget - DMS production - is 
poorly constrained, due to the difficulties associated with direct 
experimental measurement under realistic in situ conditions. 

To a first approximation, marine DMS is produced by algal and 
bacterial cleavage of DMSP released directly by phytoplankton or 
via viral and grazing activity (Gabric et al., 2001; Stefels et al„ 
2007). Apart from enzymatic DMSP cleavage, the only known 
production pathways for DMS are the enzymatic reduction of 
dimethylsulfoxide (DMSO; Spiese et al., 2009) and oxidation of 
DMSP by the hydroxyl radical (D.J. Kieber, unpublished results). 
The latter two processes are proposed to be the main source of 
DMS in marine algae that contain DMSP but do not contain 
measureable levels of DMSP lyase activity. However, it is not 
known how important these processes are as a source of dis- 
solved DMS in the water column. 

Production of DMS and DMSP by algae has various biological 
functions. DMSP is proposed to regulate osmotic pressure in the 
cells (Vairavamurthy et al., 1985) and to aid in cold tolerance 
(Kirst et al., 1991); DMSP and DMS may also act as infochemicals 
(Seymour et al., 2010) and be a part of an antioxidant mechanism 
(Sunda et al., 2002), or an overflow mechanism for excess 
reducing power (Stefels, 2000). The production of DMS by 
bacteria can be influenced by their need to assimilate sulfur from 
DMSP, which can divert sulfur away from DMS production 
(Howard et al., 2006; Simo et al., 2009). DMS loss processes 
include biological consumption by bacteria, photochemical oxida- 
tion, and removal from the ocean surface by ventilation to the 
atmosphere (Kiene and Bates, 1990; Kieber et al., 1996; del Valle 
et al., 2009). While robust approaches have been developed to 
measure microbial and photochemical removal of DMS from the 
water column, direct measurement of the biological production of 
DMS under in situ conditions is difficult because of the multiple 
ways by which DMS can be produced, including intracellular, 
extracellular, viral, and grazing-related mechanisms, and also the 
fact that the natural light field can affect DMS production (Slezak 
et al„ 2007; Gall et al., 2011). 

Toole and Siegel (2004) proposed two regimes of the pelagic 
ecosystem for DMS cycling; a stress regime and a bloom regime. 
The bloom regime refers to highly productive ocean regions, 
where the highest DMS concentrations are observed during a 
peak in primary production and chlorophyll a (Chi), which is often 
used as a proxy for phytoplankton biomass. The stress regime 
describes conditions when DMS and Chi are decoupled, and DMS 
production occurs in response to oxidative stress, as proposed by 
Sunda et al. (2002). Such conditions exist, for example, in the 
subtropical gyres, where low nutrient content and high solar 
radiation exert considerable physiological stress on marine organ- 
isms within the mixed layer. Bailey et al. (2008) conducted a 
diagnostic DMS study in the western part of the North Atlantic 
subtropical gyre (the Sargasso Sea), corresponding to a stress DMS 
regime. In contrast, the present study addressed DMS production 
in a highly productive coastal environment of the Southern 
Ocean, which we hypothesize represents a bloom DMS regime, 
where DMS production is closely correlated to Chi and DMSP 
concentrations. 

The goal of the present study was to provide a quantitative 
estimate of gross biological DMS production in the coastal 
environment west of the Antarctic Peninsula using in situ esti- 
mates of chemical and biological DMS losses and the observed 
variability in DMS concentrations and related tracers. Gross 


biological DMS production is a fundamental quantity in the 
marine DMS budget that must be well constrained in order to 
improve our understanding of DMS-climate feedbacks and their 
representation in future climate modeling efforts. The difficulties 
associated with experimental measurement of DMS production 
motivated the need for indirect estimation of gross production 
flux from more easily measured quantities. We combined water- 
column measurements collected during the DMS field campaign 
and a physical model of ocean mixing into a diagnostic calculation 
of DMS production, which relied on the difference between 
the observed and modeled DMS concentrations to infer the gross 
biological DMS production rate in the water column. 

2. Methods 

2.1. Study area 

The study area was located in the coastal waters of the 
Southern Ocean, to the west of the Antarctic Peninsula, off Anvers 
Island, near Palmer Station (64.77°S, 64.055°W). We focused on 
two stations located near the inlet to Palmer Station, B (64.78°S, 
64.075 W) and E (64.82°S, 64.04 l 'W), which are routinely 
sampled as part of the Long Term Ecological Research (LTER) 
program (Fig. 1). The LTER program has continuously operated a 
sampling network in the vicinity of Palmer Station for nearly 20 
years, starting in 1990 (http://pal.lternet.edu), providing valuable 
ancillary data to our study. The total water column depth is 
~80 m at station B and ~172 m at station E. These near-shore 
stations are affected by melt-water and runoff from the glaciers 
on Anvers Island and other small islands in the vicinity (Dierssen 
et al., 2002). Station B is located closer to the coast and the glacial 
melt-water source. Historical LTER data, available online, indicate 
that station B is characterized by lower water temperatures and 
salinities, higher Chi concentrations, and greater nutrient draw- 
down during the austral summer, compared with station E. Both 
stations experience warming and freshening of the water column 
during the austral summer but seasonal warming at station B is 
less pronounced than at station E, likely due to the glacial melt- 
water influx. A signal of deep intrusion of high-salinity oceanic 
water is evident at both stations, but is more pronounced at the 
deeper station E. 

2.2. Field measurements 

The DMS field campaign was coordinated with the LTER 2005- 
2006 summer sampling program and took place between Novem- 
ber 2005 and February 2006. Routine LTER sampling was conducted 
from a zodiac boat, approximately eveiy 2-3 d, weather and sea-ice 
conditions permitting. Conductivity-temperature-depth (CTD) 
measurements were collected with a SeaBird SBE 19 SEACAT 
instrument and processed with the SeaBird software to produce 
temperature (T) and salinity (S) profiles at 1-m resolution. Under- 
water irradiance profiles were measured with a BSI PRR-800 
Profiling Reflectance Radiometer at 16 wavelengths (Chang and 
Dickey, 2004). Discrete water samples were analyzed for Chi by 
fluorometry (Holm-Hansen et al., 1965) and for nitrate (N) by 
colorimetric spectroscopy (Knap et al., 1994). From November to 
February, in addition to routine LTER sampling, water samples 
were collected from approximately 6 depths for DMS and DMSP 
analyses. Dissolved DMS and total DMSP (DMSPt) concentrations 
were determined with a Shimadzu GC-14 gas chromatograph 
according to the procedure described in Kiene and Service (1991 ). 
In January and February, microbial consumption rate constants 
for DMS were determined with 35 S tracer methods (Kiene and 
Linn, 2000; del Valle et al., 2009). Incubations were done in the 
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Fig. 1. Map of the study area showing stations B (64.78 S, 64.075 W) and E (64.82 
LTER Palmer website http://pal.lternet.edu/sci-research/sampling-grid/). 

dark and thus our method does not take into account the effects 
of light on microbial DMS consumption, i.e., mainly inhibition due 
to UV radiation (Toole et al., 2006). Absorption coefficients of 
chromophoric dissolved organic matter (CDOM) and apparent 
quantum yields for DMS photolysis were measured in order to 
determine photolysis rate constants (Toole et al., 2003; Bailey 
et al., 2008). CDOM absorption coefficients were determined for 
the wavelength range from 210 nrn to 750 nm using a Hewlett 
Packard 8453 UV-vis photodiode array spectrophotometer. 
Wavelength-dependent quantum yields were determined 
employing a polychromatic approach outlined in Toole et al. 
(2003), using 0.2 pm-filtered sea water samples, collected on 24 
January 2006. Photochemical DMS loss was quantified using the 
35 S tracer method (Toole et al., 2004), and the UV-B (290-320 nm) 
and UV-A (320-400 nm) photon fluxes were determined by 
chemical actinometry (Kieber et al., 2007). Photolysis results 
were scaled to 1-nM DMS concentration to produce spectral 
apparent quantum yields, i.e„ moles of DMS lost per mole of 
photons absorbed by CDOM in a 1-nM DMS solution (Toole et al., 
2003). 

For this study, we focused on the period from 1 January to 
1 March of 2006, when sampling frequency was the highest due 
to favorable weather conditions and all DMS-related measure- 
ments needed to constrain the DMS budget were performed 
concurrently. The field measurements were linearly interpolated 
to a regular time and depth grid for all the following analyses. 


2.3. Diagnostic calculations 

To estimate gross biological DMS production flux from the 
available field measurements, we performed a diagnostic model- 
ing exercise based on the Reynolds-averaged statement of tracer 
conservation in turbulent flow (Holton, 1992). Conceptually, the 
average tracer concentration at any location in the flow is controlled 
by physical processes of advection and turbulent diffusion (mixing), 
and by sources and sinks of biogeochemical origin. Therefore, 
a three-dimensional conservation equation for any tracer can be 
formulated as follows: 


S, 64.041 W), which are part of the LTER basin sampling network (map adapted from 


where C is the mean component of tracer concentration, c is the 
fluctuating component of tracer concentration, U, is the mean 
velocity component of the flow, u t is the fluctuating component 
of the flow, and J is an unconstrained biogeochemical term. 
A quantitative estimate of the unknown biogeochemical processes 
can be diagnosed if all terms are constrained with observations 
and theory, except for/, which is treated as the only unknown. 
The underlying motivation behind the indirect approach to 
quantify tracer fluxes from tracer concentration data is that direct 
measurements of oceanic fluxes are difficult to make and, as a 
consequence, flux data are prone to large measurement errors. In 
comparison, tracer concentration measurements are generally 
more reliable and have a broader spatial and temporal coverage. 
In addition, the theory describing the physical processes is more 
advanced compared with the theory behind prognostic biogeo- 
chemical modeling. While the diagnostic approach is straightfor- 
ward conceptually, it becomes challenging in practical application 
because uncertainties in the diagnosed source terms are strongly 
dependent on the quality of the physical simulations and on the 
availability and quality of the data used to constrain the temporal 
and spatial tracer gradients. 

2.3.1. General model configuration 

To quantify the DMS budget at LTER stations B and E, we 
reformulated the one-dimensional form of Eq. (1) as 

— J mix 3~Jmic~l~J phot 3~Jprod< (^) 

where C [nM] is the DMS concentration, J mix [nM s _1 ] represents 
the rate of DMS gain or loss due to turbulent mixing, and the 
biogeochemical term, J, is separated into three components: the 
rate of microbial consumption, J mic [nMs -1 ]; the rate of photo- 
chemical oxidation, J phot [nM s -1 ]; and the rate of gross biological 
production, J prod [nM s 1 ], which is the objective of the diagnostic 
calculations. Atmospheric ventilation was included into the mix- 
ing term as a surface boundary condition. Horizontal advection 
was assumed negligible compared to the other terms in Eq. (2) 
(see Sections 3.1 and 3.2). 

For numerical implementation, we embedded a DMS module 
into the framework of the General Ocean Turbulence Model 
(GOTM), which is a one-dimensional physical model that, given 
the surface forcing, simulates vertical mixing in the water column 
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and solves the one-dimensional transport equations for momen- 
tum, salt, and heat energy. The model is described in Burchard 
et al. (1999) and a complete GOTM documentation is available 
online at http://gotm.net. 

The model domain depth was set to 60 m, with depth resolu- 
tion of 2 m, and the positive z coordinate was directed upward. 
The model was initialized with observed profiles of T, S, and 
DMS, and integrated from 1 January 2006 (Julian day 1 ) through 
1 March 2006 (Julian day 60) using a 30-min time step (chosen 
to ensure numerical stability of the physical simulation), with 
the biological production term excluded from the DMS budget. 
At each time step, gross biological production at a given depth was 
diagnosed from the difference between observed and modeled 
DMS concentrations: 


Nightingale et al., 2000; Blomquist et al., 2006; Vlahos and 
Monahan, 2009). The Schmidt number for DMS was calculated 
from sea-surface T according to the polynomial fit of Saltzman 
et al. (1993), and the wind speed observations were obtained 
from the automated weather system located at Palmer Station 
(Section 2.3.6). 

2.3.4. Rate of microbial DMS consumption 

The rate of microbial DMS consumption was modeled as a 
first-order process with the rate constant, k mic [s 1 ], constrained 
by in situ measurements that were linearly interpolated to match 
the model depth and time grid: 

Jmic(t.z ) = -k mic (t,z)C(t.z). (8) 


, tt-rS — C obs(f’ Z )—C mo d(U) 

Jprod\*'*Z) — ^ > w) 

where C mod [nM] is the DMS concentration calculated by the 
model, C obs [nM] is the DMS concentration produced by linearly 
interpolating measured DMS concentration profiles to the model 
vertical and temporal grid, and At [s] is the model time step. 
Following each diagnostic calculation at a given time step, model 
DMS concentrations were restored to observations at that time 
step before proceeding to the next model iteration. The model 
depth and time domains were selected based on the availability 
of observational data and excluded November and December of 
2005, when sampling was sparse due to adverse sea-ice conditions. 
The specific treatment of each rate term in Eq. (2) is discussed in 
the following sections. 


2.3.2. Rate of turbulent mixing 

The rate of vertical mixing was modeled in terms of the 
turbulent eddy diffusivity: 


JmttffA = ^t(WC) = — V(t,Z ) 


8z 


8z 


,8C(t,z) 

~8z 


(4) 


where vvc [nMms -1 ] is the vertical turbulent flux of DMS, and 
v[m 2 s'] is the vertical eddy diffusivity. In all base-run calculations, 
we used the I< profile parameterization (KPP) turbulence closure of 
Large et al. (1994) to determine the vertical eddy diffusivity. In the 
error analysis calculations described in Section 2.3.7, we used the 
k-E turbulence closure of Rodi (1987). 


2.3.3. Rate of gas exchange 

The air-sea DMS flux due to gas exchange, F gas [nMms -1 ], 
was included into the mixing term in the surface grid box and 
acted as a surface boundary condition: 


—v(t,z= 0) 


8C(t,Z = 0) 
8z 


— Fgas(t). 


(5) 


For gas flux calculations, we used the stagnant film model, in 
which, following Kettle and Andreae (2000), we assumed that 
equilibrium ocean surface DMS concentrations were negligibly 
small compared to actual surface ocean concentrations: 

F g as(t) = k w (t)C(t,z= 0), (6) 

where k w [ms -1 ] is the gas transfer velocity for DMS. To 
parameterize k w we used the formulation of Wanninkhof (1992): 

k w = yU 2 w (Sc/ 600) -05 , (7) 

where y = 0.31 (cm h -1 )/(m 2 s -2 ) is an empirical constant, U w 
[ms -1 ] is the 30-min averaged wind speed at 10m above sea 
surface, and Sc is the dimensionless Schmidt number. The 
relationship gives k„ in units of cmh -1 . In the error analysis 
calculations (see Section 2.3.7), we used a wide uncertainty 
range for k w in order to encompass the spread of other available 
gas transfer velocity parameterizations (Liss and Merlivat, 1986; 


When vertical interpolation was required, we assumed that k mic 
was equal to zero at 100 m, with the rationale that as DMS 
concentrations approached zero at that depth, microbial con- 
sumption of DMS could be assumed zero also and, consequently, 
setting k mic to zero was permissible for interpolation purposes. 


2.3.5. Rate of DMS photolysis 

Photochemical oxidation was modeled as a first-order process 
dependent on the photolysis rate constant, k ph0[ [s 1 ] : 

Jphot(Fz) = -k P hot(t,z)C(t,z). (9) 

The photolysis rate constant was determined using the algorithm 
of Toole et al. (2003): 


k P hot(.t,z)= J a CD oM(A;t)Eo(lt,z)<P(A)dl (10) 

where X [nm] is the wavelength, o C dom [m -1 ] is the CDOM 
absorption coefficient, £ 0 [Em -2 s -1 nm -1 ] is the scalar irradi- 
ance, and <P [m 3 E -1 ] is the apparent quantum yield (moles of 
DMS lost per one mole of photons absorbed by CDOM, normalized 
to DMS concentration). It should be noted that a non-SI unit, 
einstein [E], was used to represent one mole of photons. The 
integration limits were 7. 1 = 290nm and 7. 2 = 600nm. CDOM 
absorption coefficients and quantum yields were determined on 
the basis of field measurements at each station. Wavelength- 
dependent quantum yields from the measurements were fit with 
an exponential relationship, which was then used to calculate 
<F>(X) at the required wavelength resolution. Similarly, at each 
measurement time, surface wavelength-dependent o CDO m spectra 
were fit with an exponential function. The fit coefficients were 
then linearly interpolated in time to produce a C DOM(^0 fields of 
required temporal resolution. Following Toole et al. (2003), E 0 
irradiance was set to 1.2£ d (after accounting for unit conversion): 


£o = 


1.2£ d 
hvN A ’ 


( 11 ) 


where £ d [Wm -2 nm -1 ] is the measured downwellling spectral 
irradiance, h = 6.6262 x 10 -34 [J s] is Plank’s constant, v [s -1 ] is 
the frequency, and N A = 6.022 x 10 23 [mol -1 ] is Avogadro’s num- 
ber. To compute the downwelling underwater spectral irradiance 
fields at the 30-min resolution required by the diagnostic model, 
we combined the surface downwelling spectral irradiance mea- 
surements of high temporal resolution (from the spectral pyran- 
ometer at Palmer Station, see Section 2.3.6) and the spectral 
attenuation coefficients K d [m -1 ], calculated by fitting exponen- 
tial decay curves to the measured spectral underwater irradiance 
profiles of lower temporal resolution (see Section 2.2): 


E d (A,t,z) = fSE d (X, t, 0 )exp[K d (/l)z] , 


( 12 ) 


where /?=0.96 is the wavelength-independent light transmission 
across the air-sea interface, following Bailey et al. (2008). 
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Fig. 2. Summary of 30-min averaged GOTM forcing data during 1 January 2006-1 March 2006. Meteorological measurements are from the PALMOS automated weather 
system located at Palmer Station. Integrated shortwave irradiance measurements are from a spectral pyranometer located at Palmer Station. 


Table 1 

Summary of error-estimation runs and errors assigned to each input variable. 


Model run number 

Altered input variable 

Method of error assignment 

Source 

R1 

C [nM] 

Fractional error = ±10% 

Estimated from reported DMS 
measurement errors 

R2 

v [m 2 s-’] 

Used alternative turbulence closure 
scheme 

Subjective selection of Rodi (1987) 
turbulence model 

R3 

k w [ms- 1 ] 

Fractional error= ± 50% assigned to 
empirical y coefficient in the 
Wanninkhof (1992) formulation 

Subjective estimate 

R4 

kmic [s ] 

Fractional error = ±16% 

Estimated from reported 
measurement errors 

R5 

kphot [ S ] 

Fractional error = ± 50% 

Subjective estimate 


2.3.6. Model forcing 

The GOTM was forced with wind velocity, atmospheric tem- 
perature, pressure, humidity, precipitation, and downwelling 
shortwave radiation (Fig. 2). We used 30-min averaged atmo- 
spheric measurements from the PALMOS automated weather 
system located at Palmer Station and operated by the Antarctic 
Meteorological Research Center; data are available online at 
ftp://ice.ssec.wisc.edu/pub/palmer/observations. For downwelling 
shortwave radiation forcing, we used 30-min averaged integrated 
(0.285-2.8 pm) solar radiation measured with an Eppley Labora- 
tory Precision Spectral Pyranometer operated by Biospherical 
Instruments, Inc., at Palmer Station; data are available online at 
http://www.biospherical.com/nsf. 

Given that GOTM is one-dimensional, all horizontal gradients 
must be prescribed outside the model if one were to include the 
horizontal advection process into the calculation. Because the spatial 
coverage of our field measurements was limited to only two stations, 
prescribing observed horizontal tracer gradients for GOTM was not 
feasible and we chose to use a simple nudging method to account for 
the horizontal advection of T and S. The measured profiles of T and S 
were linearly interpolated to the model vertical and temporal grids 
and, at each time step, the model was relaxed towards the observa- 
tions using a one-day relaxation time-scale, i.e., the following nudging 
source term was added to both T and S conservation equations: 

Jnudg = ~ —(Ymod — Yobsi< (13) 


where r is the nudging time scale and Y represents T or S. In addition 
to horizontal advection, the nudging term potentially corrects for 
other physical processes that might play a role in the modeled 
environment but are not explicitly resolved in the physical model, 
such as a glacial melt-water source (Dierssen et al., 2002) and any 
possible tidal influence. 


2.3.7. Processing of model outputs and error estimation 

All model outputs at the original 30-min and 2-m resolution were 
averaged daily and then all days were averaged between each two 
consecutive times when the water column was sampled. Averaged 
model outputs were used in all subsequent analyses and are referred 
to as model output in the following discussion. 

To streamline analyses of the vertical structure, we divided the 
model water column into three equal-depth intervals: the 0-m to 20- 
m interval is referred to as the surface interval (SFC) and approxi- 
mately corresponds to the mixed layer; the 20-m to 40-m layer is 
referred to as the middle interval (MID); and the 40-m to 60-m layer 
is referred to as the bottom interval (BOT). The entire model water 
column, from 0 m to 60 m, is referred to as the total column (TOT). 

To provide a measure of significance of the final results, we 
selected a number of model input variables that we believed 
contributed appreciably to the final error in the model output. We 
assigned standard errors to these independent variables, combining 
available information about the measurement errors and our best 
judgment about the uncertainty in the model parameters (Table 1 ), 
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Table 2 

Propagation of errors, assigned to the input variables, into the uncertainty estimation for the dependent variables. R0 is the model base-run. R1-R5 are model error-runs as 
described in Table 1. Error notation follows Eq. (14). 


Dependent variable 

Input variables 

Partial errors in the dependent variable 

Combination of the partial errors into the total error in the 



due to each input variable 

dependent variable 

Jmix 

C 

A/ mix = Jmix,R 1 ~Jmix,R0 

(A Jmtxf = (Afabf +(4/m«) 2 


V 

A / mix = Jmix,R2 ~Jmix,R0 


Jgas 

C 

A Jgas = Jgas,R 1 ~Jgas,R0 

(AW 2 = (A/^+fAfe) 2 


kw 

AJgas = Jgas,R3 ~Jgas,R0 


Jmic 

C 

A Imic = Jmic,R 1 ~Jmic,R0 

(AW) 2 = (A/mic) 2 + (A Ct ) 2 


Is 

mic 

A/ mic Jmic,R4 ~Jmic,R0 


Jphot 

c 

^ Jphot = Jphot, R 1 ~Jphot,R0 

(Alphot) 2 = (A/phof ) 2 + < A/pCt ) 2 


kphot 

A/phof = Jphot.R 4 ~J phot.RO 


Jprod 

c 

A Jprod J prod,R\ ~J prod.R0 

(AW) 2 = (Afprod) 2 +<AW) 2 +(Aj% od f +(Afe) 2 - +<a/ p Z) 2 


V 

A Iprod = Jprod,R2~Jprod,R0 



kw 

A/ prod = Jprod,R3~Jprod,R0 



kmic 

A / prod = J prod, R4 ~J prod, RO 



kphot 

A 1 prod = Jprod,R4 ~J prod.R0 



and used a standard Taylor series expansion approach to propagate 
errors into the final results (Squires, 2001 ). For a general case of some 
variable Z that is dependent on several input variables, Z = Z(A,B,...), 
the standard error in the dependent variable, A Z, can be approxi- 
mated from the known standard errors in the input variables, AA, A B, 
etc.: : 

(AZ) 2 =(AZ a ) 2 + (AZ b ) 2 + •••, (14) 

where AZ a = (dZ/dA)AA is the change in the dependent variable 
when the input variable A is changed by the amount of its standard 
error, holding the other input variables constant, and so on. The 
implicit assumption in Eq. (14) is that the error-contributing inde- 
pendent variables are uncorrelated so that their co-variances are 
negligible (Squires, 2001 ). We calculated the model output errors due 
to each error-contributing independent variable as the difference 
between the base-run ( R0 ) and five error-runs (Rt-R5), in which the 
independent variables were altered, one at a time, as described in 
Table 1. Increasing or decreasing the input variables by the assigned 
amount resulted in approximately symmetric error in the output, that 
is, errors of approximately equal magnitude and opposite sign. A 
summary of the error propagation calculations is given in Table 2. 

3. Results 

3.1. Physical and biogeochemical characteristics 

Hydrographic observations for the modeling period are sum- 
marized in Fig. 3, with comparison between the two stations 
emphasized in the right-hand column. In Fig. 3 and all the 
following figures, the plotted error bars show one standard error 
of the means and are intended to allow visual comparison 
between the two stations. The mixed layer depth (MLD) was 
calculated as the depth at which the in situ a t exceeded the sea- 
surface a t value by 0.125 kgm~ 3 . The MLD criterion was selected 
subjectively from a range of tested criteria because it best repre- 
sented the extent of surface mixing that was visually noticeable in 
the measured hydrographic and biogeochemical tracer fields. 


The most pronounced feature emerging from this comparison was 
greater seasonal warming of the upper water column observed at 
station E, compared with station B. On average, the surface layer 
(0-20 m) at station B was slightly fresher compared with that at 
station E. These differences pointed to the influence of glacial melt- 
water at station B. Averaged over the total water column, the MLD 
was ~10m at station B and varied between 10 m and 30 m at 
station E. Hydrographic differences between the stations were 
mostly confined to the upper 20 m of the water column. 

Biogeochemical characteristics are summarized in Fig. 4, with 
comparison between the two stations shown in the right-hand 
column. For depth extrapolation, we assumed 0 nM DMS and 
DMSPt concentrations, 0 mg m~ 3 Chi concentrations, and 30 pM 
nitrate (N) concentrations at 100-m depth, which produced 
the best match with the deepest available profiles. Observed Chi 
concentrations were higher at station B than at station E at all 
depth intervals, with the most pronounced difference observed in 
the upper 20 m. Station B also had greater N drawdown at the 
surface; however, as expected, N concentrations remained rela- 
tively high even during the phytoplankton bloom. Based on the 
estimation of phytoplankton taxa using the CHEMTAX method 
(Kozlowski et al., 2011), phytoplankton species composition 
at both stations was mostly dominated by diatoms, with the 
exception of December, when cryptophytes dominated. The 
dominance of these two groups in the summer is consistent with 
the analysis of the 13-yr time series (1995-2007) on the shelf 
along the western Antarctic Peninsula (Kozlowski et al., 2011). 
DMSPt concentrations were clearly lower at station B compared 
with station E in the upper 20 m and were similar below 20 m 
depth. Averaged over the entire 60 m water column, DMSPt con- 
centrations were 50 + 8 nM at station B and 74 + 10 nM at station E 
(Table 3). Despite the observed differences in DMSPt, Chi, and to a 
lesser extent N, the DMS distributions were very similar at the two 
stations, with total-column averages of 3.3 ± 0.6 nM at station B and 
3.2 ± 0.8 nM at station B. As emphasized in the right-hand column of 
Fig. 4, the two stations were virtually identical in terms of the average 
DMS concentrations. This observation suggested that it was reason- 
able to assume that horizontal DMS gradients in the region were 
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Fig. 3. Temperature, salinity, and density at sea-level pressure measured at stations B and E during 1 January 2006-1 March 2006. Interpolated fields are shown in color 
and observation times are indicated by light vertical lines to illustrate the extent of interpolation. Calculated mixed layer depth is shown as a thick white line (see Section 3.1). 
The last column compares the physical characteristics by depth interval (TOT: 0-60 m, SFC: 0-20 m, MID: 20-40 m, BOT: 40-60 m), averaged over the entire modeling period. 
Error bars give one standard error of the means estimated as follows: the time series were averaged by depth intervals defined above; the mean and the standard deviation were 
calculated for each depth-averaged time series; the standard error of the mean was estimated as the standard deviation divided by the square root of the number of elements in 
the time series. (For interpretation of the references to color in this figure legend, the reader is referred to the web version of this article.) 


close to zero, which then can serve as a justification for neglecting the 
horizontal advection term in the DMS budget, as we have done. 

3.2. Diagnosed DMS budgets 

The time and depth-interval averaged DMS budget parameters for 
the base-run, i.e„ measurement-based rate constants for microbial 
consumption ( k mic ) and photolysis (k ph ot ), and model-based vertical 
eddy diffusivity (v) and gas transfer coefficient (k w ), are summarized 
in Fig. 5. Similar to the station characteristics discussed above, 
differences in the budget parameters between stations B and E were 
observed mostly in the surface layer (0-20 m). On average, k mic was 
slightly greater at station B. Average k pho t was non-zero only in 
the surface layer for both stations, and was clearly greater at 
station E, which, according to our analysis, had lower light 
attenuation. Dierssen et al. (2002) reported that turbidity and 
radiance reflectance (ratio of upwelling radiance to downwelling 
irradiance) were about two times higher at station B compared 
with station E, based on the analysis of the LTER data from 1991 
through 1999. The authors suggested that highly scattering 
minerogenic particles released by melt-water cause the water to 
be optically distinct from typical conditions and appear to sink 
out rapidly offshore, leading to dramatically different optical 
conditions at stations B and E. Average model v and k m were 
practically indistinguishable at the two stations. 

A summary of the diagnosed budgets for stations B and E is 
presented in Fig. 6. In the surface layer, the plotted rate of mixing 
represents a combination of two processes: turbulent mixing and 
gas exchange. At both stations, ventilation was about 8% of the 
total losses in the surface interval. Mixing can lead to both gains 


and losses of DMS, depending on the direction of eddy movement 
and the DMS concentration gradient. When averaged, mixing 
rates were very close to zero in the interior and acted as a very 
minor loss term for DMS in the surface interval, accounting for 
less than 1% of the total losses at both stations. Microbial 
consumption was the dominant loss term in both station budgets, 
in agreement with previous field experiments in Antarctic waters 
(del Valle et al., 2009) and in other oceanographic settings (Kieber 
et al., 1996; Vila-Costa et al., 2008). DMS loss to microbial con- 
sumption was stronger at station B, while photolysis loss was 
stronger at station E. In the surface interval, microbial consumption 
accounted for 80% and 66% of the total losses at stations B and E, 
respectively, while photolysis represented 12% and 26% at stations B 
and E, respectively. The DMS production rate was approximately 
equal at both stations. It is to be noted that, since microbial DMS 
consumption rates were measured in dark incubations, and this 
process has proved susceptible to UV radiation-mediated inhibition 
(Toole et al., 2006), the biological DMS loss may have been over- 
estimated, leading to overestimation of the diagnosed gross biolo- 
gical DMS production rates. 

From our diagnostic analysis, the scale for all of the source 
terms in the DMS budget was on the order of 1-10 nM d -1 . We 
used this scale information to provide additional justification for 
omitting the horizontal advection term in the DMS budget. The 
horizontal advection rate,y ad „, in Eq. (1) can be scaled as 

Jad,~U X' (15) 

where U is the horizontal velocity scale and (AC/I) is the scale for 
the horizontal gradient in the DMS concentration. According to 
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Fig. 4. DMS, DMSPt, Chi, and N concentrations measured at stations B and E during 1 January 2006-1 March 2006 period. Interpolated fields are shown in color and the 
observations are indicated in white to illustrate the extent of interpolation. The last column compares the biogeochemical characteristics by depth interval (TOT: 0-60 m, 
SFC: 0-20 m, MID: 20-40 m, BOT: 40-60 m), averaged over the entire modeling period. Error bars give one standard error of the means as described in the Fig. 3 caption. 
(For interpretation of the references to color in this figure legend, the reader is referred to the web version of this article.) 


Table 3 

Diagnosed rates of DMS production, normalized production rates, and concentrations of DMS, DMSPt, and Chi from the present study of the coastal Southern Ocean (B and E) a in 
comparison with results from the Sargasso Sea study (Cl and A2) by Bailey et al. (2008). 


Interval 

[m] 

Station 

DMS 

[nM] 

Jprod 

[nM d -1 ] 

DMSPt 

[nM] 

Jprod/DMSPt 

[(nM d _1 )/nM] 

Chl 

[mg m~ 3 ) 

Jprod/Chl 

[(nM d _1 )/(mg m -3 )] 


B 

3.3 ± 0.6 

3.1 ±0.6 

50 ±8 

0.06 ± 0.02 

6 ± 1 

0.5 ±0.1 

0-60 

E 

3.2 ± 0.8 

2.7 ± 0.6 

74 ±10 

0.04 ± 0.01 

3.8 ± 0.7 

0.7 ± 0.2 

TOT 

Cl 


0.73 ± 0.09 

15.5 

0.047 ± 0.006 

0.063 

12 ± 1 


A2 


0.9 ±0.12 

10.4 

0.087 ±0.014 

0.044 

20 ±3 


B 

4.7 ± 0.8 

5 ± 1 

88 ±11 

0.06 ± 0.01 

9 ± 2 

0.5 ±0.1 

0-20 

E 

4.6 ±0.1 

5 ± 1 

137 ± 19 

0.04 ± 0.01 

5.7 ± 0.9 

0.9 ± 0.3 

SFC 

Cl 


0.64 ± 0.09 

10.8 

0.058 ± 0.008 

0.040 

15 ± 2 


A2 


0.78 ± 0.09 

7.5 

0.104 ±0.012 

0.033 

22 ±3 


B 

3.1 ±0.6 

2.7 ± 0.5 

39 ±9 

0.07 ± 0.02 

6 ± 1 

0.5 ±0.1 

20-40 

E 

3.1 ±0.8 

1.7 ±0.4 

53 ±10 

0.03 ± 0.01 

3.4 ± 0.8 

0.5 ± 0.2 

MID 

Cl 


0.91 ±0.14 

14.5 

0.062 ±0.010 

0.055 

15 ± 3 


A2 


1.14 ± 0.16 

9.5 

0.120 ±0.017 

0.040 

28 ±4 


B 

1.9 ±0.4 

1.4 ±0.3 

24 ±7 

0.06 ± 0.02 

4.0 ± 0.9 

0.4 ± 0.1 

40-60 

E 

1.8 ±0.5 

1.0 ±0.2 

32 ±8 

0.03 ± 0.01 

2.3 ± 0.6 

0.5 ± 0.2 

BOT 

Cl 


0.64 ±0.10 

21.3 

0.030 ± 0.005 

0.095 

7 ± 1 


A2 


0.78 ±0.15 

14.1 

0.055 ±0.011 

0.059 

13 ± 3 










a For normalized rates, the error bars give the combined error due to the propagated uncertainty in J pr0 d (determined as described in Section 2.3.7) and the standard 
error of the mean DMSPt and Chi (determined as described in Fig. 3 caption). The squared combined fractional error was estimated as the sum of squared fractional errors 
in the components, assuming that the error-contributing variables were uncorrelated. 
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Fig. 5. Summary of DMS budget parameters used in the base-run budget calculations for stations B and E. Measured rate constants for microbial consumption (k mic ) and 
photolysis [k phot ), and modeled vertical eddy diffusivity (v) were averaged by depth interval (TOT: 0-60 m, SFC: 0-20 m, MID: 20-40 m, BOT: 40-60 m); modeled gas 
transfer coefficients ( k w ) are for the surface model grid box only. All parameters were averaged over the entire modeling period (1 January 2006-1 March 2006). Error bars 
give one standard error of the means determined as described in the Fig. 3 caption. 
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Fig. 6. Diagnosed DMS budgets for stations B and E. The DMS production rate 
( Iprod), combined rate of mixing and gas exchange (J, n j X ), microbial consumption 
rate (/ mlc ), and photolysis rate (J phot ) are averaged by depth interval (TOT: 0-60 m, 
SFC: 0-20 m, MID: 20-40 m, BOT: 40-60 m) and over the entire modeling period 
(1 January 2006-1 March 2006). Error bars give estimated uncertainty in model 
outputs calculated by propagation of errors from the model input variables 
(see Section 2.3.7). 

Smith et al. (1999) and Klinck et al. (2004), surface currents in the 
region are ~0.01 m s -1 . We used the average DMS concentration 
difference between stations B and E from our analysis (Fig. 5) to 
approximate AC~0.1 nM, and the horizontal distance between 
the stations, L, is approximately 5000 m. We estimated the scale 
for the horizontal advection rate at J od „ ~ 0.01 nM d _1 , which is 
between two and three orders of magnitude smaller than the 
other source terms in the DMS budget. 

3.3. Diagnosed rates of DMS production 

Total column-averaged time series of DMS production were 
compared with the biogeochemical parameters in Fig. 7. The 




patterns emerging from this comparison were similar for both 
stations. Notwithstanding the uncertainty, there were two peaks 
in the production rates - around Julian days 23 and 40 - reaching 
close to 7 nM d _1 , while background rates were less than 5 nM d _1 . 
Peaks in DMS production were reasonably well aligned with 
peaks in DMSPt and Chi and with minima in N concentrations. 
A strong storm event occurred around January 23, with wind 
speeds reaching up to 25 ms -1 and the atmospheric pressure 
decreasing by 22 hPa, down to 973 hPa, over the course of the day 
(Fig. 2). The wind was predominantly from the north, changing 
from northeasterly to northwesterly as the storm progressed. 
There were clear increases in observed DMS and DMSPt concen- 
trations (Fig. 4) and a spike in diagnosed DMS production (Fig. 7) 
during the storm, suggesting a possible storm-induced stress 
related DMS production. Peaks in DMSPt concentrations were 
slightly higher at station E, which is located further offshore and 
could have been affected more by the storm because of less 
sheltering by the surrounding land masses. However, because the 
error bars in diagnosed production rates overlapped, we were 
unable to detect a difference in the DMS production at the two 
stations in response to the storm. Another possible cause for the 
observed increase in DMS and DMSPt concentrations could have 
been upwelling of water from near the sediment-water interface, 
where sediment-derived sources of DMS and DMSP might have 
existed, as has been shown, for example, for the coastal zones of 
Northern Europe (Nedwell et al., 1994; Belviso et al., 2006) and 
for the Ross Sea, Antarctica (DiTullio et al., 2000). As emphasized 
in Fig. 4, our coverage of the deeper layer of the water column 
was sparse with few measurements collected below 30 m depth 
because the study focused on the upper water column processes; 
thus, we cannot conclusively attest whether the sediment source 
of DMS and DMSP was present in our study area. Because our 
stations are fairly shallow, extending to the depths of ~80m 
(station B) and ~172 m (station E), it is possible that a sediment 
source may affect upper water column concentrations during 
a deep mixing event. However, such transport of bottom water 
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Fig. 7. Time series of the diagnosed DMS production rate (Jprod ), DMSPt, Chi, and N for stations B and E between 1 January 2006 and 1 March 2006. Time series are averaged 
over the entire model water column (0-60 m). For J prod , error bars give estimated uncertainty in model output calculated by propagation of errors in the model input 
variables (see Section 2.3.7). For DMSPt, Chi, and nitrate (N) concentrations, error bars give one standard error of the means estimated as the standard deviation divided by 
the square root of the number of elements in each vertical profile. 
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Fig. 8. Diagnosed DMS production rates (] pr od ), production rates normalized to DMSPt (J prod /DMSPt), and production rates normalized to Chi (J pr0 d/Chl) for stations B and E. 
Rates were averaged by depth interval (TOT: 0-60 m, SFC: 0-20 m, MID: 20-40 m, BOT: 40-60 m) and over the entire modeling period (1 January 2006-1 March 2006). For 
J P r 0ll , error bars give the estimated uncertainty in model output calculated by propagation of errors from the model input variables. For normalized rates, error bars give the 
combined error due to the propagated uncertainty in J pro{t (see Section 2.3.7) and the standard error of the mean DMSPt and Chi (determined as described in Fig. 3 caption). 
The squared combined fractional error was estimated as the sum of squared fractional errors in the components, assuming that the error-contributing variables were 
uncorrelated. 


would have required overturning of the seasonal stratification 
and mixing of the entire water column, which was unlikely during 
stable summer stratification. In addition, based on the available 
deep profiles, both DMS and DMSPt concentrations showed a 
clear decrease with depth (Fig. 4) and the depth of the mixed 
layer did not exceed 30 m at either station (Fig. 3), which also 
implies that a sediment source was unlikely. 

A depth-interval summary of diagnosed production rates, and 
production rates normalized to DMSPt or Chi concentrations is 
given in Fig. 8. Averaged over the entire water column, diagnosed 
DMS production rates at stations B and E were 3.1 + 0.6 and 
2.7 + 0.6 nM d _1 , respectively. At both stations, the average 
production rate of DMS in the surface mixed layer was about 
two times higher than in the interior. Averaged over the entire 
water column, the ratio of DMS production to Chi at station E 
was slightly higher than at station B (0.7 + 0.2 and 0.5 + 0.1 
(nM DMS d _1 )/(mg Chi m~ 3 ), respectively). On the other hand, 
DMS production rates normalized to DMSPt standing stock 
were slightly higher at station B at all depths. Averaged over 


the entire column, the DMSPt-normalized DMS production was 
0.06 + 0.02 (nM DMS d _1 )/(nM DMSP) at station B and 0.04 + 0.01 
(nM DMS d- 1 )/(nM DMSP) at station E. 

Averaging both stations, the DMS production rate in the water 
column that we inferred through our diagnostic calculation was 
2.9 + 0.4 nM d -1 . Combining this estimate with the average observed 
DMS concentration of ~3 nM, we estimated the DMS replacement 
time - DMS concentration divided by gross production rate - to be on 
the order of 1 d, which was similar to the average DMS replacement 
time estimated by Gall and Simo (2010) in productive Arctic waters 
in the summer, and at the lower end of turnover times estimated for 
microbial DMS consumption reported in a wide range of oceano- 
graphic settings (Simo and Pedros-Alio, 1999). 


4. Discussion 

Observed peaks in the diagnosed production rates at stations B 
and E approximately followed the maxima in Chi concentrations 
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(Fig. 7), consistent with the bloom regime of DMS production 
concept proposed by Toole and Siegel (2004). Moreover, we found 
a significant positive correlation between observed Chi concen- 
trations and the diagnosed DMS production rates (Fig. 9), indicat- 
ing that phytoplankton biomass is a strong driver of DMS 
production in the Antarctic Peninsula study area. 

We compare our results with the results of other DMS studies 
conducted in the Southern Ocean region and with a diagnostic 
DMS study in the contrasting biogeochemical environment in the 
North Atlantic subtropical gyre by Bailey et al. (2008). 

4.2. Comparison with previous studies in the Southern Ocean 

Berresheim et al. (1998) reported DMS concentrations in the 
surface waters near Palmer Station ranging between 0.7 and 
3.7 nM during January-February 1994, which was slightly below 
the values we observed in January-March 2006 (5.6 + 3.4 nM at 
station B and 5.5 + 4.2 nM at station E). In a different location in 
the Southern Ocean, a similar temporal evolution of DMS, DMSPt, 
and Chi concentrations to our January-March 2006 observations 
was described by Gambaro et al. (2004) for the Gerlache Inlet 
(Terra Nova Bay, Ross Sea, Antarctica), with two clear maxima. 
Our observed DMS concentrations, however, were approximately 
an order of magnitude lower than in the Gerlache Inlet study. In 
another study in the vicinity of Gerlache Inlet, Kiene et al. (2007) 
reported surface water measurements along three transects 
between 49° S and 76° S, in the New Zealand sector of the Southern 
Ocean, from November 2003, December 2004, and November 
2005. In open waters, DMS concentrations ranged between 
0.6 and 3.2 nM, but reached up to 30 nM at the northern 
boundaries of the seasonal sea-ice (63° S-68° S) and in the 
northern Ross Sea. Good agreement with our observations was 
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Fig. 9. Scatter plots of diagnosed DMS production rates (J prod ) and observed Chi 
concentrations averaged over the surface interval (0-20 m) for stations B and E. 
Also shown are the correlation coefficients and the associated p-values. 


also found in the study by Jones et al. (2010) in the south-west 
Indian Ocean sector of East Antarctica, from 30° E to 80° E 
and south of 62° S, from January to March 2006. Surface DMS 
concentrations ranged from non-detectable to 48 nM, with a 
mean value of 10 nM, which was slightly higher than our average 
concentrations. Further east, also in the Australian sector of the 
Southern Ocean (63° E-162° E, 40° S-69° S), Curran and Jones 
(2000) reported a monthly average surface water DMS concentra- 
tion of 7.9 nM for the spring and summer seasons from 1991 to 
1995 in the seasonal ice zone, which is comparable to average 
DMS concentrations observed in our study. 

DMSPt concentrations reported by Kiene et al. (2007) ranged 
from 12 to 52 nM in the open waters and reached up to 95 nM at 
the seasonal sea-ice boundary and in the Ross Sea. While our 
DMSPt concentrations agree on average, on several occasions we 
detected much higher DMSPt concentrations, reaching up to 
400 nM at station E. This may reflect much higher phytoplankton 
biomass near Palmer Station (up to 30pgC/i/L _1 ) compared to 
the Ross Sea (~7 pg Chi L _1 ), but it may also reflect potential 
losses of DMSPt during preservation of Ross Sea samples contain- 
ing colonial Phaeocystis antarctica (del Valle et al., 2010). Lower 
mean surface DMSPt concentrations, ranging from 10 to 52 nM, 
with a mean of 29 nM, were also reported for the Australian 
sector by Jones et al. (2010). 

With respect to DMS losses, our estimates of biological 
DMS consumption in the surface layer (0-20 m) ranged from 
2 to 7nMd~' (4.2 + 0.8) at station B and from 1 to 8nMd~' 
(3.6 + 0.7) at station E. These rates are in the same range as those 
measured in the Ross Sea, up to 8.8 nM d _1 (Kiene et al., 2007; del 
Valle et al., 2009). Curran and Jones (2000) reported an average 
gas flux from the seasonal ice zone of 49 pmol m~ 2 d -1 , for an 
average wind speed of 11.7 m s _1 . In our study, the average wind 
speed was 3.5 m s _1 , and the estimated gas fluxes were slightly 
lower, ranging between 2 and 32pmolm~ 2 d~ 1 (Fig. 10). In 
addition, Curran and Jones (2000) estimated the DMS emissions 
in the Antarctic region of the Southern Ocean (south of the Polar 
Front, about 55° S), excluding contributions form sea-ice, at 
85 Gmol S yr ’, representing 10% of the mean global emission 
estimate of 880 Gmol S yr -1 provided by Lana et al. (2011) from 
6% of the global ocean surface area. These results suggest that the 
Southern Ocean and, in particular, the seasonal ice zone, is an 
important DMS source, emphasizing the need of continued study 
of DMS dynamics in this region. 

In summary, the DMS concentrations observed in our study are 
in good agreement with other studies in the Southern Ocean. The 
same is true about modeled air-sea DMS flux and microbial DMS 
consumption rates. Observed DMSPt concentrations appear higher 
than what is reported in other studies. To the best of our knowledge, 
there are no published estimates of gross biological DMS production 
in the Southern Ocean and, therefore, we are unable to compare our 
diagnosed rates with observations from literature. 
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Fig. 10. Time series of the modeled DMS air-sea gas flux for stations B and E during 1 January 2006-1 March 2006. A positive flux is from the ocean to the atmosphere. 
Error bars give the estimated uncertainty in model output calculated by propagation of errors in the model input variables (see Section 2.3.7). 
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4.2. Comparison with a previous study in the North Atlantic 
Subtropical gyre 

We compared our results with the results of a diagnostic DMS 
study in the North Atlantic subtropical gyre by Bailey et al. (2008), 
referred to as B08 hereafter (Table 3). The B08 study used a 
similar modeling approach to estimate gross biological DMS 
production in two eddies in the Sargasso Sea, a downwelling 
ocean region characterized by low nutrient concentrations and 
extremely low primary productivity in the summer. Therefore, 
B08 and the present study of the coastal Southern Ocean region 
describe two contrasting biogeochemical ocean environments. 
Average gross biological DMS production in the water column 
estimated in the present study (2.9 + 0.4 nMd -1 ) was about 
three times greater than in B08. The DMS replacement time 
estimated in the present study was on the order of 1 d, which 
was shorter than in B08, where the reported replacement time 
was between 4 and 9 d, indicating faster DMS turnover in 
the Antarctic coastal waters despite much lower temperatures 
(1-2 °C at Palmer Station vs. ~30°C in the Sargasso Sea). Gross 
production in B08 was the greatest below the mixed layer, 
whereas in the present study gross production was the greatest 
within the mixed layer. Compared with our estimates, gross 
biological DMS production normalized to Chi was about an order 
of magnitude higher in B08. This likely reflected the differences in 
the composition of phytoplankton assemblage and the low 
phytoplankton biomass in the Sargasso Sea surface waters in 
summer. Average DMSPt:CM ratios were an order of magnitude 
higher in B08 ( ~ 100 nmol pg -1 ) than in the present study 
(~10nmol pg -1 ), again due to low chlorophyll biomass, as the 
Antarctic DMSPt concentrations were 5-fold higher, on average, 
than in B08, and up to 400-fold higher in surface waters. We 
found, however, a striking similarity in production rates normal- 
ized to DMSPt concentrations, suggesting a strong underlying 
similarity between two contrasting biogeochemical environments 
in terms of DMSPt-DMS conversion. The significance of this 
finding of near constancy of DMSPt-DMS conversion ratio across 
pelagic biomes is that it could provide a means to use DMSPt 
measurements to extrapolate gross DMS production, which is 
extremely difficult to measure experimentally. 

B08 conducted a comparison of their DMS production esti- 
mates with those reported in literature for other ocean regions 
and pointed out that the range in reported DMS gross production 
rates spans nearly three orders of magnitude, from ~0.6 to 
500 nM d \ After normalizing to DMSPt, however, the range covers 
only about an order of magnitude, from 0.05 to 0.7 (nM DMS d -1 )/ 
(nM DMSP). Averaging the results of B08 with the results of the 
present study, our estimate of mean gross biological DMS production 
normalized to DMSPt concentration, independent of the biogeo- 
chemical environment, is 0.06 ± 0.01 (nM DMS d -1 )/(nM DMSP). 
This represents a conversion of about 6% of the standing stock of 
DMSP into DMS per day, a value consistent with the bulk of DMSP- 
sulfur being diverted to non-DMS metabolites (Kiene and Linn, 2000; 
Howard et al., 2006). 

4.3. Inference about global marine DMS production 

Lastly, we used our estimate of average DMSPt-normalized 
DMS production to make a first, albeit crude, observation-based 
estimate of the magnitude of global marine gross biological DMS 
production. Kettle et al. (1999) estimated global-median surface- 
ocean dissolved and particulate DMSP concentrations at 9.7 nM 
and 20.5 nM, respectively, from which we inferred, by summa- 
tion, a global-median surface-ocean DMSPt concentration of 
30 nM. Reported dissolved and particulate DMSP medians were 
computed from data binned onto a 1 -degree grid by calendar 


month (about 1000 values for each quantity). We chose to use the 
median rather than the mean as a measure of the central 
tendency of the global DMSPt concentration in the calculation 
because Kettle et al. (1999) emphasized that the distribution of 
the binned data was strongly non-Gaussian, skewed towards 
lower values. We subjectively estimated that the median DMSPt 
concentration has an error of ± 50%. Multiplying this estimate by 
our DMSPt-normalized DMS production gives a global average 
surface ocean gross DMS production estimate of 2 + 1 nMd -1 . 
Assuming that this DMS production occurs in the upper 100 m 
and using the global ocean area of 3.6 x 10 s km 2 , we arrive at a 
global ocean gross DMS production rate of 24 ± 13 Tmol S yr '. 
Using the Lana et al. (201 1 ) global estimate of the sea-to-air DMS 
flux, 0.88 Tmol S yr -1 (range from 0.55 to 1.08 Tmol S yr -1 ), DMS 
gas emissions represent 4 + 2% of marine gross DMS production 
globally. In a global modeling study, Kloster et al. (2006) simu- 
lated DMS gas emissions of 0.88 Tmol S yr -1 , in agreement 
with Lana et al. (2011), and marine DMS production of 
11 Tmol S yr -1 , which is about half of our estimate. Despite these 
differences, these independent estimates both arrive at the same 
conclusion that, at the global scale, the DMS flux from the ocean 
to the atmosphere is only a small fraction of the marine DMS 
production. 


5. Summary and conclusions 

We estimated gross biological DMS production in the coastal 
environment west of the Antarctic Peninsula using in situ esti- 
mates of chemical and biological loss terms in the DMS budget 
and the observed variability in concentrations of DMS and related 
tracers. The calculations were based on water-column mea- 
surements collected during the DMS field campaign and a 
physical model of ocean mixing. Averaged over the entire water 
column, diagnosed gross biological DMS production rates were 
3.1 + 0.6 nM d -1 at station B and 2.7 ± 0.6 nM d -1 at station E. 
Production rates were at a maximum in the mixed layer at both 
stations. The ratio of DMS production to Chi was 0.7 + 0.2 and 
0.5 + 0.1 (nM DMS d -1 )/(mg Chi m 3 ) at stations B and E, respec- 
tively. The DMSPt-normalized gross biological DMS production 
was 0.06 ± 0.02 (nM DMS d -1 )/(nM DMSP) at station B and 
0.04 + 0.01 (nM DMS d -1 )/(nM DMSP) at station E. It should be 
emphasized that our field campaign did not include measure- 
ments of turbulent mixing and gas exchange processes; therefore, 
these processes add considerable uncertainty to our results 
because their treatment was entirely model-based in our calcula- 
tions. Future studies should include direct measurements of 
mixing and gas exchange processes in order to have data that 
could be used for model calibration. 

A comparison of our results with the results of a diagnostic 
DMS study in the North Atlantic subtropical gyre by Bailey 
et al. (2008) revealed that the production rates normalized to 
DMSPt concentrations were very similar in the two contrasting 
biogeochemical ocean environments. The significance of this 
finding of near constancy of DMSPt-DMS conversion ratio 
across pelagic biomes is that it could provide a means to use 
DMSPt measurements to extrapolate gross biological DMS pro- 
duction, which is difficult to measure in situ. We estimated 
that mean gross biological DMS production normalized to 
DMSPt concentration, independent of the biogeochemical envir- 
onment, is 0.06 ± 0.01 (nM DMS d -1 )/(nM DMSP). Combining 
this estimate with currently available values for global-median 
surface-ocean DMSP concentrations (Kettle et al., 1999) suggests 
a global ocean gross biological DMS production rate of 
24 ±13 Tmol Syr -1 . 
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